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We provide an analytical description ol the dynamics ol an atom in an optical lattice using the 
method of perturbative adiabatic expansion. A precise understanding of the lattice-atom interac- 
tion is essential to taking full advantage of the promising applications that optical lattices offer in 
the field of atom interferometry. One such application is the implementation of Large Momentum 
Transfer (LMT) beam splitters that can potentially provide multiple order of magnitude increases 
in momentum space separations over current technology. We also propose interferometer geometries 
where optical lattices are used as waveguides for the atoms throughout the duration of the inter- 
ferometer sequence. Such a technique could simultaneously provide a multiple order of magnitude 
increase in sensitivity and a multiple order of magnitude decrease in interferometer size for many 
applications as compared to current state-of-the-art atom interferometers. 
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I. INTRODUCTION 



Atom interferometry has opened new frontiers in precision metrology. Highly sensitive gravimeters, gravity gra- 
diometers, and gyroscopes have been constructed, and promising work has been done to integrate these sensors into 
a robust apparatus that can operate outside the laboratory with applications in inertial navigation and geodesy 
[iHS]. Moreover, atom interferometers have been used to make competitive measurements of the fine structure con- 
stant [101 lll| . Since atom interferometric measurements of the fine structure constant do not assume the validity 
of Quantum Electrodynamics (QED), while determinations of the fine structure constant based on measurements 
of the electron magnetic moment do make this assumption, comparison between the results of these two methods 
provides a stringent test of QED [T2|, [T3] . In addition, an experiment to test Einstein's Equivalence Principle with 
unprecedented precision is underway |14) . and atom interferometric gravitational wave detectors offer the possibility 
to study gravitational radiation in frequency ranges complementary to LIGO and LISA |15| . 

Atom interferometers have traditionally relied on matter gratings or light pulses to act as beam splitters and mirrors 
for matter waves, with atomic wave packets traveling freely between these interaction zones. Light-pulse schemes using 
either Raman pulses (where the internal state of the atom is changed) or Bragg pulses (where the internal state of the 
atom remains unchanged) have been implemented, such as those described in |16H19) . For a number of applications 
of light-pulse atom interferometers, such as measurements of gravity and rotation, the sensitivity is proportional to 
the separation in momentum that can be attained between the two arms |20) . In measurements of the fine structure 
constant, the sensitivity scales as the square of this separation [TU]. Therefore, significant efforts have been devoted 
to the development of Large Momentum Transfer (LMT) beam splitters. LMT beam splitters achieving momentum 
splittings of 24:hk using multi-photon Bragg pulses have recently been demonstrated [17]. However, the required laser 
intensities to make significant improvements on this result may prove to be prohibitive |21j . In contrast, LMT beam 
splitters that use several two-photon Bragg pulses or a multi-photon Bragg pulse of relatively small order to separate 
the two arms of the interferometer in momentum space, followed by the acceleration of one of the arms with an 
optical lattice, could potentially provide multiple order of magnitude increases in attainable momentum separations 
with relatively modest laser intensity requirements. 

An atom interferometer that uses this method has been successfully operated in a proof of principle experiment 
(with a maximum demonstrated momentum splitting of 12hk) |22) . In a separate experiment, an atom interferometer 
with lOhk LMT beam splitters has been realized using a similar technique [23] . Alternatively, both arms of the 
interferometer could be simultaneously accelerated in opposite directions by two different optical lattices after the 
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FIG. 1: (color online) Schematic of a guided lattice interferometer. The trajectories of the two arms are guided by two separate 
optical lattices for the duration of the interferometer sequence, resulting in a fully confined interferometer. The two paths 
shown in the figure are the trajectories of the two portions of the atomic wavepacket, which respectively correspond to the two 
arms of the interferometer. 

initial splitting. Using this second scheme, an interferometer with 24/ifc LMT beam splitters that achieves 15% contrast 
and an individual beam splitter that provides an 88hk momentum separation have been demonstrated p4j . 

The utility of atom interferometry hinges upon the ability to precisely calculate the phase accumulated along the 
different arms of an interferometer |25H27| , of which the phase acquired during interactions of the atoms with light is 
an important component. Indeed, the phase obtained by an atom during a Raman or Bragg pulse is well-understood 
[T1|2T1|28]. Analogously, in order to take full advantage of the potential of lattice beam splitters, we must have a 
detailed understanding of the phase evolution of an atom in an optical lattice. In this paper, we provide a rigorous 
analytical treatment of this problem. To our knowledge, such a treatment has not been previously presented in the 
literature. 

Based on this analysis, we propose atom interferometer geometries in which optical lattices are used to continuously 
guide the atoms, so that the atomic trajectories are precisely controlled for the duration of the interferometer sequence, 
with a different lattice guiding each arm of the interferometer (as illustrated in Fig. [T]). We point out here a distinction 
in terminology between a lattice waveguide and a lattice beam splitter. Here, a lattice waveguide is the use of 
a lattice to continuously control the trajectory of an arm of an atom interferometer. We note that two separate 
lattice waveguides can independently control the two arms of an interferometer, or a single lattice waveguide can 
simultaneously control both arms. In contrast, a lattice beam splitter is an interaction of relatively short time (in 
comparison to a waveguide) with the primary purpose of splitting the arms of the interferometer in momentum space 
rather than providing continuous trajectory control. The underlying physics behind lattice waveguides and lattice 
beam splitters is the same, and they can be treated with a common formalism. A single lattice waveguide that 
simultaneously transfers IQOOhk of momentum to the two arms of a Ramsey-Borde interferometer has been previously 
achieved in [llj . However, to our knowledge, our idea of using optical lattice waveguides to create a fully confined 
atom interferometer has not been previously considered. 

Our analysis indicates that these lattice interferometers will offer unprecedented sensitivities for a wide variety of 
applications and that they will be able to operate effectively over distance scales previously considered too small to 
be studied by precision atom interferometry. For example, one particularly interesting configuration involves using 
two optical lattice waveguides to continuously pull the two arms of the interferometer apart, subsequently holding the 
two arms a fixed distance from each other in a single lattice waveguide that is common to the two arms, and then 
using two lattice waveguides to recombine the arms. Such a configuration could be used, for instance, as a gravimeter. 
The sensitivity of lattice interferometers is illustrated by the fact that, given the experimental parameters stated in 
[T] (10^ atoms/shot and 10~^ shots/s), a shot noise limited lattice gravimeter whose arms are separated by 1 m for 
an interrogation time of 10 s has a sensitivity of 10~^^ g/Hz^/^. We perform phase shift calculations for these lattice 
interferometers using the theoretical groundwork formulated in this paper, and we discuss how lattice interferometers 
can both exceed the performance of conventional atom interferometers in many standard applications and expand the 
types of measurements that can effectively be carried out using atom interferometry. 

The paper is organized as follows. Sec. II. describes the Hamiltonian for an atom in an optical lattice in the 
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different frames we use in the paper. Sec. III. discusses the phase evolution of an atom in an optical lattice under 
the adiabatic approximation. Sec. IV. introduces the formalism of perturbative adiabatic expansion to calculate 
corrections to the adiabatic approximation, and Sec. V. applies this formalism to calculate phase corrections to a 
lattice beam splitter. Sec. VI. proposes a number of interferometer geometries that make use of lattice manipulations 
of the atoms. The main results of the paper are Eqs. (27) and (28), which show how to obtain analytical corrections 
to the lowest order phase shift estimates. These corrections are surprisingly large, and understanding them is vital to 
realizing the full accuracy of the sensor geometries proposed in Sec. VI., as well as other geometries utilizing optical 
lattice manipulations of the atoms. For example, the gravitational wave detector proposed in |15j will likely make 
use of lattice beam splitters and/or waveguides. Previously, the phase evolution induced by lattice manipulations 
was not sufRciently well-understood to allow for a detailed design of the atom optics system or an estimation of the 
corresponding systematic effects. 



II. THE HAMILTONIAN IN DIFFERENT FRAMES 



An optical lattice is a periodic potential formed by the superposition of two counter-propagating laser beams. 
Atoms can be loaded into the ground state of the lattice by ramping up the lattice depth adiabatically, and the lattice 
can then be used to impart momentum to the atoms and/or to control the atoms' trajectories. Optical lattices are 
thus a useful tool for atom optics. 

We begin our discussion of the lattice-atom interaction by finding a useful form for the Hamiltonian. As is typical 
for many applications of atom interferometry, to minimize decoherence we assume that we work with atomic gases 
dilute enough so that the effects of atom-atom interactions are negligible. We first consider the Hamiltonian in the lab 
frame, where for now we assume a vertical configuration with constant gravitational acceleration g so that we have a 
gravitational potential given by mgx. We expose the atom to a superposition of an upward propagating beam with 
phase </>up(i) and a downward propagating beam with phase ^down(i), which couples an internal ground state \g) to 
an internal excited state |e). The two-photon Rabi frequency is = ^up(*)_^°wn(*) ^ .^^j^gj-g .^^g jg^ riup(f) denote the 
single-photon Rabi frequency of the upward propagating beam, 57down(i) denote the single-photon Rabi frequency of 
the downward propagating beam, and A denote the detuning from the excited state. We depict the physical setup 
in Fig. [2] Making the rotating wave approximation and adiabatically eliminating the excited state as is standard 
procedure pO] . we obtain the following Hamiltonian-where the periodic term in the potential arises from a spatially 
varying AC stark shift and where k is the magnitude of the wave vector of the laser beams [29l |30j : 



i?Lab = ^ + "ih^it) sin^ 
2m 



mgx 



(1) 



Note that where the difference between the frequency of the upward propagating beam and the frequency of the 
downward propagating beam is denoted by Auj{t), we will have the relation A(j){t) = (f>up{t) — 4>downit) = Jq Auj{t')dt' + 

0up(O)--(/'down(O). For a given A0(i), the lattice standing wave will be translated by DLab(0 = ^tk ^ in the x direction 
from the origin. Thus, the velocity of the lattice in the lab frame is: 



-Lab(t) = |^Lab(<) = ^ (2) 



and we rewrite the lab frame Hamiltonian as: 



i?Lab - ^ + ^hilit) sin^ [kx - fcDLab(i)] + mgx (3) 
2 m 

In order to most readily describe the dynamics of an atom in an accelerating optical lattice, it is useful to work 
in momentum space. The mgx term that appears in the lab frame Hamiltonian makes such an approach difficult, 
especially when considering non-adiabatic corrections to the phase shift. However, we can change frames by performing 
a unitary transformation in order to obtain a Hamiltonian that is easier to handle analytically. In the end, we will 
see that approaching the problem from the point of view of dressed states provides a convenient Hamiltonian for 
our purposes. We consider the transformation procedure from the lab frame to the dressed state frame in Appendix 
A, where we also introduce an intermediate frame that freely falls with gravity (which we call the freely falling 
frame) . We note that the general form of the unitary transformations considered in Appendix A as well as the specific 
transformations to the different frames we consider can also be found in the Appendix of 30J. 
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FIG. 2: (color online) The physical setup for applying a periodic potential. We expose atoms to counter-propagating laser 
beams with respective single-photon Rabi frequencies Q,^p and r^down. The lasers are detuned from the transition between the 
atom's internal ground state and excited state so that the atom's external momentum states are coupled through two-photon 
transitions, creating an effective lattice potential. 

It is convenient to absorb the initial velocity wo of the atom in the lab frame into the dressed state frame, so that 
velocity vq in the lab frame corresponds to velocity zero in the dressed state frame. The Hamiltonian in the dressed 
state frame is, as derived in Appendix A: 



HTys = ^- (wLab(i) +gt- vo)p + 2hn{t) sin2 (kx) (4) 

Now, we will show how working in momentum space allows us to represent i?DS as an infinite dimensional, discrete 
matrix. This matrix is discrete because the optical lattice potential term, Vosin^(fci), only couples a momentum 
eigenstate \p) to the eigenstates \p + 2hk) and \p — 2hk) [30]. For the moment, we will examine the evolution of 
individual eigenstates of the dressed state Hamiltonian iJ^s • These eigenstates reduce to single momentum eigenstates 
\p) when = 0. The knowledge of how each of these eigenstates evolves under i^DS will allow us to describe the 
dynamics of an entire wavepacket. For the moment, we will only consider momentum eigenstates corresponding to an 
integer multiple of 2hk, since we have boosted away the initial velocity vq of the atom in the lab frame. We note that 
it is always possible to transform to a particular dressed state frame in which a given momentum eigenstate in the 
lab frame corresponds to zero momentum in that dressed state frame. The results we derive here can thus be readily 
generalized to arbitrary momentum eigenstates in a wavepacket, as we discuss in greater detail in Appendix B. 

We consider a discrete Hilbert space spanned by the momentum eigenstates \2nhk) for integers n, so that we can 
express any vector in this Hilbert space as: 



\^{t))= ^ Cn{t)\2nhk) (5) 

n— — oo 

Since this Hilbert space is discrete, it is natural to adopt the normalization convention that {2mhk\2nhk) = 5mn- 
When considered as an operator acting on this discrete Hilbert space, i?DS can be written as pTl 15^ : 



^discrete ^ }__L \2nhk) {2nhk\ - J2 {vLe.hit)+gt-vo)i2nhk)\2nhk){2nhk\ 

n— — OO n— — oo 

°° n(f) 

- {\2nhk) {2{n ~ l)hk\ + \2nhk) {2{n + l)hk\) 



(6) 



where we drop the common light shift. Now, it is convenient to introduce the recoil frequency u!r = ^ = |^ and the 
recoil velocity «,. = —. In order to make our notation as compact as possible, we will be interested in the quantity 
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a{t) = "Lab (tHg* ^'o ^ which is the velocity of the lattice in the dressed state frame in units of w^. We can express the 
second term of ifg'l""*'' in a useful way by noting that (i;Lab(i) + gt - vo)2nhk = ^^-^(*)+9*-''o y^2nhk = Ana{t)Er. 
Furthermore, we define Cl{t) = ^.We can now write the discrete Hamiltonian in a simplified form: 

oo 

^discrete ^ ^2nnk) {2nhk\ - na{t) \2nhk) {2nhk\ - n{t) {\2nhk) {2{n - \)hk\ + \2nhk) (2(n + l);ifc|)] (7) 

11 — — OO 

The matrix elements of this Hamiltonian are: 

77™„ = {2mhk\ iJ^I"-*^ \2nhk) = AE^ [(n^ - na{t)) (5„„ ~ n{t) (<5„,„+i + (5™,„_i)] (8) 

In matrix notation, the Schrodinger equation for the discrete Hilbert space takes the form ih^^{t) = H'^{t), where 

we let H be the matrix whose element in the mth row and nth column is given by Hmn and we let \E'(i) be the column 
vector whose nth entry is c„(t). 

III. PHASE EVOLUTION UNDER THE ADIABATIC APPROXIMATION 

Now that we have determined the Hamiltonian matrix H, we have the appropriate machinery in place to describe 
the phase evolution of an atom in an optical lattice. We consider the process in which momentum is transferred to the 
atom through Bloch oscillations. Reference [3D] provides a thorough and insightful description of Bloch oscillations 
in a number of different pictures. Given our choice of Hamiltonian, we work in the dressed state picture, which is 
discussed in Sec. IV. B. of [30], making extensive use of Bloch's theorem, the concept of Brillouin zones, and the 
band structure of the lattice [21]. As in the previous discussion, we consider the evolution of single eigenstates of the 
dressed state Hamiltonian, noting that we can easily generalize our results to the case of a wave packet of finite width, 
as we address in Appendix B. 

Initially, we consider the system to be in a momentum eigenstate. First, we adiabatically ramp up the lattice 
depth by increasing the laser power so that we load the system into an eigenstate of the Hamiltonian. For the 
purposes considered here, we want the system to enter the ground eigenstate (corresponding to the zeroth band of 
the lattice). In order for this to be achieved, a resonance condition must be met, which states that the velocity of 
the lattice must match the velocity of the atom to within Vr- The loading will be adiabatic if the adiabatic condition 

{1\H\0) <C j^" is satisfied, where |0) and |1) respectively denote the ground state and the first excited state 

of the Hamiltonian [22] [30] . This condition will be easier to meet near the center of the band (where the velocity of 
the lattice is identical to the velocity of the atom), because the energy gap ei — Eq between the zeroth band (which 
corresponds to ground state) and the first band (which corresponds to the first excited state) becomes smaller as the 
velocity difference between the lattice and the atom becomes larger (which corresponds to moving toward the border 
of the first Brillouin zone). The resonance condition is discussed further in Appendix C. 

In the lab frame, the atom accelerates under gravity, increasing the deviation between its velocity and the lattice 
velocity during the loading process. In the freely falling and dressed state frames, in which gravity is boosted away, 
this corresponds to the lattice accelerating upward while the atom remains at rest. This effect can negatively impact 
the loading efficiency if the loading sequence is sufficiently long so that the accrued velocity difference becomes a 
significant fraction of Vr ■ In such a scenario, the effect can be ameliorated by accelerating the lattice in the lab frame 
to fall with the atom, which corresponds to the lattice velocity remaining constant in the freely falling and dressed 
state frames. Furthermore, we note that in the case of a lattice LMT beam splitter, the lattice should only be resonant 
with one arm of the interferometer, so that negligible population from the other arm is affected. Otherwise, the signal 
could be distorted by multi-path interference, causing a systematic error in the estimation of the interferometer phase 
shift. Conditions for when the negative effects of off- resonant lattices can be avoided can be estimated using the 



Hamiltonian matrix for an off-resonant lattice given in Eq. ( 33 1 . We discuss off- resonant lattices quantitatively and 
in more detail in Sec. VI. 

After the adiabatic loading of the atom into the ground state of the Hamiltonian, the frequency difference between 
the laser beams is swept to accelerate the lattice, periodically imparting momentum to the atom in units of 2Hk. 
In the dressed state frame, this phenomenon can be understood in terms of avoided line crossings, which occur 
because the coupling between the atom and the laser beams lifts the degeneracy at the crossing points. We refer the 
reader to Fig. 10 of [HU] for a clear illustration of these avoided crossings. As the frequency difference is swept so 
that the system passes through the avoided crossings, the system remains in the ground state of the dressed state 
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FIG. 3: (color online) The left panel shows a numerical simulation of a lattice acceleration in momentum space that transfers 
lOfifc of momentum. The lattice depth and velocity are shown as a function of time in the right panel, where in this particular 
case the relevant parameters are Tramp = 120{4ijJr)^^ , Taccci = 16(4aj,.)~^, fimax = and Omax = 10 (corresponding to a 

final lattice velocity of 10 Vr)- First, we adiabatically ramp up the lattice to a depth of 19.5Er so that the lattice is loaded into 
the ground state of the dressed state Hamiltonian. Subsequently, we accelerate and ramp down the lattice, leaving the atom 
in a single momentum eigenstate. 

Hamiltonian as long as the process is adiabatic. Consequently, at each of the avoided crossings, the momentum of 
the atom increases by 2tik, which corresponds to a Bloch oscillation. Finally, after the acceleration of the lattice, the 
frequency difference is held constant while the lattice depth is adiabatically ramped down, delivering the system into 
the momentum eigenstate \pi + 2Nhk) , where pi is the momentum before the Bloch oscillations and N is the number 
of Bloch oscillations. 

Fig. [3] depicts the lattice depth and velocity as functions of time for the process described above and shows a 
numerical simulation of a particular instance of this process: the adiabatic loading of the lattice from an initial state 
\OHk), the transfer of lOhk of momentum through 5 Bloch oscillations, and the ramping down of the lattice to deliver 
the system into the final state \10hk). 

Since under the adiabatic approximation we assume that the atom always stays in the ground state of the Hamil- 
tonian, the phase (/)(i) of the atom evolves as follows: 



eo{t')dt' 



(9) 



where £o{t) is the instantaneous ground state eigenvalue of the Hamiltonian and ^0 is the initial phase of the atom. 
In addition to Eq. (|9|, there is also a Berry's phase term [33]. However, this term is zero for a linear external 
potential. Therefore, there is no contribution from the Berry's phase under the semiclassical approximation, as long 
as the external potential is treated as linear-we discuss the validity and ramifications of this approximation at the 
end of Appendix A. We note that any such contribution would arise from the residual external potential terms of 
the dressed state Hamiltonian that are non-linear in x (we collectively denote these terms as V in Appendix A and 
explain why they can often be neglected). 

We now consider how the eigenvectors and eigenvalues of H change with a, which we recall is the dimensionless 
velocity of the atom in the dressed state frame. Say that the eigenvalues of H(a^fl) are given by e„(Q;,r2) with 
corresponding eigenvectors ^'„(q!, fi), where the index n runs from to cxd. We choose to index the eigenvalues so 
that e„(a,ll) denotes the the nth eigenvalue labeled in order of increasing value. Moreover, we let c'f\a,n) be 

The transformation 



the jth element of the column vector ^„(a, 57), so that ^„(q!, fl)) = J2i= 



,c)->{a,n) \2jhk). 



properties of the eigenvectors and eigenvalues under changes in a can be deduced from Bloch's theorem. It can be 
shown that when the lattice velocity is increased by 2vr = , which corresponds to a being increased by two, while 
J7 is kept fixed, the new eigenvectors can be obtained through the following relation [3T] : 



{a - 



2,n) = c'-"],{a,n) 



(10) 
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This simply represents a shift of the wavefunction in momentum space by 2hk, which is exactly what we expect, since 
increasing the lattice velocity by 2vr corresponds to undergoing a single Bloch oscillation. The dependence of the 
eigenvalues on the lattice velocity can be expressed as follows: 



(11) 



where p„(a, fi) is periodic in a such that Pn{a + 2m, fl) — Pri(a, il) for integer m holds for all a and fl) vanishes 
when the condition a — 2m holds. This periodicity follows directly from Bloch's theorem, since increasing a by 2m 
corresponds to increasing the lattice velocity by 2mvr, meaning that a and a + 2m are at the same point in the 
first Brillouin zone. Note that the dependence of the eigenvalue on can be calculated using the truncated matrix 
approximation discussed in Sec. IV. The relevance of this dependence to the phase shift of an interferometer and how 
this dependence varies with momentum are discussed in Sec. V. and Appendix B. 

The first term in Eq. ( 11 ) has a simple physical interpretation. Where f Lattice (^) is the velocity of the lattice in the 
dressed state frame (so that fLattice(i) = ctVr), note that: 



Era^ 



h k WLattice(i) 



1 



2m 



hk 



mvi^atticcity 



(12) 



which is simply the kinetic energy of an atom traveling along a classical trajectory defined by the motion of the lattice. 
The e„(0,r2) and p„(q!, f2) terms represent the band structure of the lattice, with the Pn{a,Cl) term accounting for 
the bands deviating from being flat. For interferometer geometries in which lattices act as waveguides for the atoms, 
the net contributions to the phase shift from the eo(0,r2) and pQ(a,fl) terms in the ground state energy and from 
corrections to the adiabatic approximation are often negligible, as explained in the following sections. In this case, 
the phase difference between the two arms of the interferometer is given by (assuming that the two arms arrive at the 
same endpoint): 



„arinl 



{t)dt 



„arm2 



{t)dt 



"^Latti. 



-'Lattice 



(13) 



where T is the time elapsed during the interferometer sequence. 

Observe that this expression for the phase difference can be obtained by assuming that the lattice potential acts as 
a constraint that forces the atoms in each arm to traverse the classical path traveled by the lattice guiding that arm. 
In this case the phase shift is just the difference of the respective action integrals over the two classical paths, as we 
would expect from the Feynman path integral formulation of quantum mechanics 134! . Since the lattice is the only 
potential in the freely falling and dressed state frames, the action integrals yield Eq. ( jlS] ) (for an insightful treatment 
of the applications of path integrals in atom interferometry, we refer the reader to 35 ). The terms that we have 
neglected in Eq. ( 13 ) embody corrections to the simple picture of the lattice as a force of constraint arising from the 



quantum nature of the motion (e.g., a small portion of the population leaving the ground state of the lattice), which 
can sometimes be important. However, our simple picture provides physical intuition into the lattice phase shift and 
is often sufficient to derive quantitative results. 

To summarize, the eigenvalues of the lattice consist of a kinetic energy term, a term that depends only on the 
lattice depth, and a term that is periodic in a, and the eigenvectors transform under a simple shift operation when a 
is changed by 2m for integer m. These properties are a direct result of Bloch's theorem. The symmetries that we have 
discussed allow us to conclude that if, for a given fl, we know the eigenvalues and eigenvectors of the Hamiltonian for 
all a within any range [Q;o,ao + 2], we can subsequently determine the eigenvalues and eigenvectors for arbitrary a. 
This result will prove to be useful from a computational standpoint, since the dynamics of the system are completely 
described by the solution within a finite range of a. 



IV. CALCULATING CORRECTIONS TO THE ADIABATIC APPROXIMATION USING THE 
METHOD OF PERTURBATIVE ADIABATIC EXPANSION 



We now present the method of perturbative adiabatic expansion |21j to determine corrections of arbitrary order 
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to the adiabatic approximation. We note that the particular adiabatic approximation that we correct here refers to 
the adiabatic evolution of the ground state of the dressed state Hamiltonian, rather than the adiabatic ehmination 
of the excited state during the Raman process, which is treated in [3T]. The corrections we consider will always be 
present to some extent, since lattice depth and velocity ramps occurring over a finite time can never be perfectly 
adiabatic. In addition, non-adiabatic corrections can be caused by perturbations arising from laser frequency noise 
and amplitude noise. Although our analytical and numerical computations indicate that the contribution of non- 
adiabatic corrections to the overall phase shift will be highly suppressed for many interferometer geometries, it is 
important to have a generalized framework with which to treat these corrections in order to determine when they are 
important and to precisely calculate them when necessary. 

During the course of this derivation, it is convenient to parameterize the Hamiltonian, eigenvalues, and eigenvectors 
using the single variable t rather than the two variables a{t) and (l{t). Note that much of our discussion will follow a 
similar outline as the proof of the adiabatic theorem in |33] . A more detailed version of the derivation presented here 
can be found in |31j . 

For all times t, we can express any state vector ^'(i) in Hilbert space as a linear combination of the instantaneous 
eigenvectors ^'ri(t) of the dressed state Hamiltonian, where in general the coefficients of each eigenvector can vary in 
time. The instantaneous eigenvectors satisfy the relation i/(t)^'„(t) = e„(i)vl'„(t). Choosing coefhcients with a phase 
iPn{t) = — ^ en{t')dt' factored out, we can write: 



(14) 



To simplify matters further, we choose the phase of ^^(i) so that each element of ^n(i) is real for all t and varies 
continuously with t, which we can do because the particular Hamiltonian matrix H(t) we consider is a real-valued, 
Hermitian matrix. 

When applied to Eq. (14 1, the Schrodinger equation gives us the relation: 



oo 



b^(t) = - ^6„(i)*](t)^^e'[v."(*)--^Wl (15) 

n=0 

We note that the inner products (t) ^"^qI^^ are zero for the case of a linear external potential, as verified numerically, 
and we will thus drop them from the sum. These inner products are closely related to the Berry's phase, because 



integrating ^] {t) ^^gt*^ with respect to time gives the Berry's phase for the jth eigenvector. Eq. (15) provides 
us with a relation which we could directly use to perform adiabatic expansion. But to see more clearly how the 
adiabatic expansion series relates to the rate at which the Hamiltonian changes in time, we will express bj{t) in a 

more transparent form. As long as e„(<) — £j{t) ^ 0, we can express ^] {t) ^'^g/*"' in terms of a matrix element of H{t) 
and an energy difference |31) : 



dt 



We can therefore write Eq. ( 15 ) as: 



£n{t)-£j{t) 



(16) 




K{m{t) 



dt 



^ m{t)H(mjt) ., ,,,, 



(17) 



where Spjt) is the set of all n such that n ^ j and £„(i) = ej{t) and SNoit) is the set of all n such that e„(i) ^ £j{t)- 



Eq. il7h illuminates the rationale behind the adiabatic approximation. Under the adiabatic approximation, we 



assume that H(t) and hence also its eigenvectors vary slowly enough in time so that the conditions ^](i)i?(i)^'„(<) 

{6„{t)-ej{t)f ^^^^ ^ ^^^^^^ vl/t(t)^^^L(!) < 1 (for e„(t) = e.j{t) and where At is the time scale of the 

approximation) hold. The righthand side of Eq. (171 can therefore be approximated as zero. Then, all the coefhcients 
bj{t) are constant in time. 
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To compute higher order corrections, we employ the method of adiabatic expansion, which mathematically follows 
in the spirit of the Born approximation. To zeroth order, we take the coefficients bj{t) to be constant as dictated by the 
adiabatic approximation. To obtain the first order co rrec tions to these coefficients, we substitute the constant zeroth 
order coefficients b'"^^ into the righthand side of Eq. ( 15 ) with n ^ j and integrate to find the first order coefficients 

h^p{t). We can repeat this process recursively, substituting the coefficients b^l\t) into Eq. (15) to calculate the 

coefficients bj {t) and continuing until we know the coefficients to the necessary precision. This method provides a 
way to construct a series expansion for each hj (t) . We adopt a matrix notation for the terms in this series to facilitate 
the discussion. The first order correction to bj{t) consists of contributions from each nonzero bn "^ where n ^ j, and 
we depict the contribution from bn ^ as C„_>.j. So to first order, we can write: 



bf\t)^bf^+J2c^^, (18) 

where: 



The second order solution for bj{t) will then be 



where: 



C™^„^, ^ -^*C™^„^t(i.)5Mle'K(0-^.(0]di' (21) 

under the implicit assumption that the time variable associated with Cm-^n is appropriate for the context in which it 
appears. The calculation of corrections of higher order is discussed in Appendix D. 

To find the eigenvectors and eigenvalues that we need to calculate the terms that make a non-negligible contribution 
to the expansion, we must approximate the infinite dimensional Hamiltonian matrix as a finite dimensional truncated 
matrix. At the end of Sec. III., we concluded that the problem of determining the eigenvalues and eigenvectors for 
all a reduces to finding the eigenvalues and eigenvectors for a range a € [Q;o,ao + 2] for arbitrary a^. In addition, 

it suffices to calculate the inner products {t) ^'^g/*"' just in this range of a, which follows from the symmetry 

^j(a + 2m, Vl) = 'i/j(a, fj) ^'^"^^'^) fQj- integer m [Slj. To make the calculation less cumbersome, we 

can look at the range a G [—1, 1]. For $7 not too large and a in this range, the eigenvectors with lower energies are 
populated almost entirely by momentum eigenstates \2mfik) with relatively small |m|. This is the case because for 
a in the range, the diagonal elements of the Hamiltonian will be smallest for values of m close to zero. We note 
that for = 0, the diagonal elements are the eigenvalues. In the limit of 11 — > 0, each eigenvector will consist of 
only a single momentum eigenstate, where in general eigenvectors corresponding to momentum eigenstates with m 
closer to zero will have lower eigenvalues. Increasing will allow the lower eigenvectors to spread out in momentum 
space to a certain extent, but this will not change the fact that the lower eigenvectors will be linear combinations 
of momentum eigenstates corresponding to smaller values of |to|. As discussed previously, the eigenvectors we care 
about for calculational purposes will be those with eigenvalues closer to the ground state eigenvalue. We can thus 
consider a truncated (2n + 1) x (2n + 1) Hamiltonian matrix centered around m = 0, where we choose n to be large 
enough so that for the particular dynamics being described, a sufficient number of eigenvectors and eigenvalues can 
be calculated. 



V. AN EXAMPLE OF PERTURBATIVE ADIABATIC EXPANSION: CALCULATING THE 
NON-ADIABATIC CORRECTION TO THE PHASE SHIFT EVOLVED DURING A LATTICE BEAM 

SPLITTER 

We illustrate the above method by calculating phase corrections to a lattice beam splitter. In this example, we 
consider the case where two optical lattices of the same depth but different accelerations are used to separate the two 
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arms of the interferometer (after an initial momentum space splitting is achieved through Bragg diffraction) . We note 
that the analysis here is equally applicable to the situation where two separate optical lattice waveguides are used to 
address the two arms of the interferometer, an example of which is illustrated in Fig. [7j To calculate the phase shift 
for applications in precision measurement, we need to determine the non-adiabatic correction to the phase difference 
between the two arms that accrues during the beam splitter. In practice, we do this by first calculating corrections 
to the ground state coefficient bo{t) and then evaluating how these corrections affect the phase difference between the 
arms. We note that the dominant contribution to the phase difference will come from the zeroth order term as given 
inEq. 

For this example, we consider the situation shown in Fig. |3] , where the interaction of the atoms with the lattice is 
divided into three distinct parts. From t = to t = Tlamp we ramp up the lattice, from t = Tramp tot = T^amp + Thecal 
we accelerate the lattice, and from t = Tiamp + Taccci to t = 2Ti.amp + ^accei = Tgnai we ramp down the lattice. For the 
sake of simplicity, the ramps are chosen to be symmetric so that the lattice depth decrease ramp is the time reversed 
lattice depth increase ramp. 

We assume that initially all of the population is in the ground state, so that fej'^-' — Sqj, and we make the lattice 
depth and velocity ramps adiabatic enough so that almost all of the population remains in the ground state. In order 
to find the non-adiabatic correction to the phase shift, we determine the non-adiabatic correction to the phase of the 
ground state for each arm. 

Since to lowest order only the ground state is populated, the leading corrections to foo(^) will come at second order. 
Using the formalism in Eq. ( 20 1 , the second order correction will be: 



6f w-C^ = E^o— 



(22) 



The largest contribution comes from Co-i.i-i.o(^)i it is this term on which we focus. Because the ground state is 
non-degenerate, Eqs. (17) and (21 1 give us: 



Jo ei(ii) - £o(^i) 

\ Jo Aeio(t2) J -Aeioih) 

(23) 

where we define Mio(t) = '^\{t)H {t)'^o{t) ~ '^o{t)H {t)^ i{t) (note that the two matrix elements are equal because 
we choose the eigenvectors to be real) and A£io(i) = eo(i) — 

We examine the ultimate contribution of Co->i-i.o(21inai) to 6q (Tfinai)- Since during the ramp up and ramp down 

stages Mio{t) depends only on fl but not on d, some portions of Co->i-i.o(2finai) will be common to both arms of the 
interferometer, because we assume that the lattice interaction processes for the two arms differ only in the magnitude 
of the lattice acceleration. We denote these common terms as (?ramp- The remaining terms depend on the lattice 

acceleration and will thus differ between the arms. There will be a term gmixod that depends both on f2 and a. 
However, it can be shown that under the assumption that the lattice depth decrease ramp is the time reversed lattice 
depth increase ramp, this term is zero |31j . Finally, there will be a term gaccei that depends quadratically on a and 

is not explicitly dependent on fl. We note that (jaccci implicitly depends on the maximum lattice depth fimaxi which 
can be seen by the fact that J^max affects what value the energy gap Aeio takes on during the acceleration stage (a 
deeper lattice leads to a larger energy gap). We can thus write: 



&0 (^final) ~ 1 + (/ramp + .9accel (24) 

In calculating .gaccei, it is useful to note that Mio(i) takes on a convenient form during the acceleration stage. 
Recalling the form of the Hamiltonian matrix from Eq. (jsj), we observe that H(t) will be a diagonal ma- 
trix with matrix elements H„in{t) — —4:Erna{t)Smn- We can thus write Mio(t) = a{t)Aio{t), where Aio{t) = 
-AEr Y,n (^o(i)) is a weighted dot product. 

In order to more clearly illuminate the general points we are illustrating with this example, we make the simplifying 
assumption that d(<) is constant throughout the acceleration stage. Moreover, we assume that the lattice is deep 
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enough so that Aio{t) and Aeio{t) are also constant during the acceleration stage, which is an accurate approximation 
for typical experimental situations. During the acceleration stage, we respectively denote these constant quantities 
as a, AiQ, and Aeiq^^^^^. Note that these assumptions, along with the assumption of mirror symmetry between the 
ramp up and ramp down stages, are certainly not necessary to carry out the calculation. They only serve to make 
the final result take a particularly simple form that provides physical insight into the process. In the absence of 
these assumptions, the calculation will be only slightly more complicated and can easily be performed. We note in 
particular that the mirror symmetry assumption is not stringent, for even when this symmetry is largely violated, the 
5mixcd term is typically an order of magnitude or more smaller than the gaccci term, as we verify by estimating the 
relevant integrals in Eq. (23). If needed, (/mixed can be calculated by evaluating these integrals. In addition, we note 
that the treatment given in this example can readily be generalized to the case where the lattice depth and velocity 
are changed simultaneously. Performing the necessary integrals, we find that: 



— —ih 



(Aeio_„,)^ 



accel 



(Aeio_,)^ 



(25) 



We note that for |f/ramp + 5accei| <C 1, we can solve the problem by employing adiabatic expansion over a single time 
interval. However, there may be times when we must divide the problem into multiple parts, as discussed in Appendix 
D. For typical experimental parameters, the term proportional to Taccei in .gaccci will dominate both the second term 
in gaccci and the (7ramp term. We have verified that the (/ramp term (which embodies the non-adiabatic loading of the 
lattice) is typically much smaller than the second term in gaccoi by estimating the integrals in Eq. ( 23 ) that correspond 
to gramp and by checking these estimates numerically. Thus, we can express the condition l^ramp + <?accci| ^ 1 as: 



(A£io_,) 



2 ^acccl 



< 1 



(26) 



This condition will often hold, since in many experimentally relevant cases the acceleration time or acceleration will 
be sufficiently small. 

We now show how to determine the correction to the phase shift between the two arms arising from the non- 
adiabatic correction gramp + Saccei in the case where this correction is small. As we recall from Eq. ( 14 ), the coefficient 
of the ground state eigenvector at time Tfinai is ^o(2finai)e*'^"*^^'''"'''-' ~ (1 + .giamp + 5accci) e*'^''*^^""'''^ Now, note that the 
argument of any complex number z — |z|e''^''s(^) can be written as arg(z) = ^ In [pr] . The non-adiabatic correction 
to the phase of the coefficient of the ground state eigenvector for an arm with acceleration a is thus: 



(1 + Ifframpl e'"^s(g— p) + Iffaccell e'^^S^g— ')) 
(1 + IfframpI 6"* ^>-g(9-»p) ~t- Iffacccll e"* -'■g(ff--i) ) 

~ IfframpI sin (arg [(/j-amp]) + l5accci|sin(arg[5accoi]) 

(27) 

where we have Taylor expanded to first order in small quantities. To calculate the correction to the phase shift 
between two arms, we take the difference </)corrcction(aarmi) - </>corroction(aarm2)- The Ifframpl sin (arg [giramp]) term is 
common to both arms. The non-adiabatic correction to the phase difference between an arm with acceleration a and 
an unaccelerated arm is: 



^correction ((^) — i^ 



correction V^corrcction (0) « Isacccil sin (arg [ffaccoi]) (28) 



Eqs. (27 1 and ( [28^ provide us with a means to calculate the leading correction to the phase shift, and comparison 
with numerical results for a variety of experimentally conceivable lattice depth ramps shows excellent agreement. A 
comparison of the adiabatic expansion method with numerical calculations is illustrated in Fig. |4] In particular. Fig. 
I^jshows S(f){a) as a function of a for an acceleration time of 16(4^,.)"^ with Omax = The numerical values shown 
m the figure come from a numerical simulation of the Schrodinger equation. The non-adiabatic correction to the 
phase shift between the two arms of the interferometer for arbitrary arm acceleration differences is '/'correction (c^armi) — 
0correction('iarm2) = '^^^'(aarmi) ~ i50(<iarm2)- Notc that the leading Order results depicted in Fig. [4] will often be 
sufficient, but it may sometimes be necessary to calculate higher order corrections, examples of which can be found 
in i31j- 
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FIG. 4: (color online) Non-adiabatic correction to the phase difference between an arm that is accelerated by an experimentally 
plausible lattice beam splitter and an arm that is not accelerated (which is in our notation 5(j>{a)), as calculated using a 
simplified adiabatic expansion method in which we keep only leading terms versus numerical simulations of the Schrodinger 
equation. The curve represents the prediction made by the adiabatic expansion method, while the red dots represent the 
numerical results. The lattice depth ramps and acceleration time are identical to those of the acceleration sequence shown 
in Fig. [sj with an acceleration time of 16(4tJr)~^ and with f2max = In addition to the leading second order term, we 



keep the leading fourth order correction term in Qaccei, —ft' 



2 A^Al 



2(Ae 



calculated in (31] , which becomes significant for larger 



accelerations (e.g., for an acceleration of 6 (^) Vr, this term is smaller than the leading second order correction by a factor of 
~ 10). We neglect corrections arising from the term flramp, for these corrections are common to both arms of the interferometer 
to lowest order. Even the simple approximation used to obtain the curve agrees remarkably well with the simulations (with an 
rms deviation of 4 x 10"** radians), and we note that we could easily improve this approximation by including more terms in 
the adiabatic expansion series. As expected, we observe that the correction scales quadratically with acceleration. 



In an experimental implementation, the two arms of the beam splitter are addressed by two different lattices. 
Therefore, any imbalance in the depths of the two lattices will lead to a phase error between the arms. Where the 
intensities of the two beams forming a given lattice are la and 7^, we note that the lattice depth is proportional to 
the product -v/Tq x ^/h, meaning that intensity imbalances lead to lattice depth imbalances. We calculate that once 
a lattice is ram ped up, the eo(0, ^) term in the expression for the lowest eigenvalue of the dressed state Hamiltonian 
from Eq. (Ill is much larger than the po^a^fl) term. Thus, for an interaction lasting from time ti to time ^2, the 
dominant contribution to the phase error, which we denote as (/Jhedanca arising from the imbalance in the lattice depths 
is: 



-^balance 



£o(0i ^arml) — £o (Oj ^aim2) 



dt 



(29) 



where we note that eo(0,n) can be calculated using the truncated matrix approximation discussed in Sec. IV. In 
order to put Eq. (291 in a more convenient form for making order of magnitude estimates, we use the fact that for 
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|f2| > 0.25, eo(0,r2) ^ —8Er\fl\. This result is verified by direct comparison with the values of eo(0, 17) obtained with 
the truncated matrix approximation. Recalling that Q = ^j- , it is convenient to rephrase this statement in terms of f2 
as follows: for |fi| > 2L0r, eo(0, ^) ~ — This range of |0| is of considerable experimental interest and is amenable 
to simple approximation. However, if needed, smaller lattice depths can be treated using the general relation given 
in Eq. (29) and the truncated matrix approximation. Substituting the above result into Eq. (29), we obtain: 



^balance ^ / ( I ^^arml I " I f7arm2 I ) = ( ( I ^^arml |> - ( I ^7arin2 I ) ) (^2 " ^ 1 ) (30) 
Jti 

where (|i7armi|) and (|r2arm2|) respectively denote the average values of |riarmi| and |fiarm2| between time ti and time 
t2- Note that fluctuations in the difference between |riarmi| and |f2arm2| that occur at frequencies that are large with 
respect to the beam splitter time t2 — ti will largely average out, suppressing their net effect on the phase shift. 
Furthermore, in many geometries, arm 1 will be addressed by one pair of laser beams, that we call lattice A, during 
the splitting of the arms and then will be addressed by a second pair of laser beams, that we call lattice B, during the 
recombination of the arms. Conversely, arm 2 will be addressed by lattice B during the splitting stage and by lattice A 
during the recombination stage. In such a geometry, the effect on the phase shift of a constant offset in depth between 
lattice A and lattice B will cancel, and the effect on the phase shift of fluctuations in the depth difference between 
lattice A and lattice B that occur at low frequencies with respect to the time scale of the interferometer sequence will 
be highly suppressed. Also, as we discuss in the following section, we will often be interested in the difference in the 
phase shifts of two interferometers in a differential configuration. If both of the interferometers remain well within 
the Rayleigh ranges of the laser beams so that beam divergence is a small effect, any lattice depth imbalance will be 
largely common to the two interferometers, suppressing its effect on the phase shift difference by orders of magnitude. 



VI. APPLICATIONS: ATOM INTERFEROMETERS USING OPTICAL LATTICES AS WAVEGUIDES 

Light-pulse atom interferometer geometries have had tremendous success in performing many types of high-precision 
measurements. However, in many cases, we would like to be able to push the capabilities of atom interferometry by 
making more precise measurements using spatially compact interferometers. Atom interferometers that use optical 
lattices as waveguides for the atoms offer the potential to make such measurements attainable. In such a scheme, we 
can use an initial beam splitter composed of multiple Bragg pulses, a multi-photon Bragg pulse, or a hybrid Bragg 
pulse/lattice acceleration scheme as described in [22-,-,24j to split the arms of the interferometer in momentum space. 
We can then control each arm independently with an optical lattice. We will once again use Bragg pulses during 
the TT-pulse and final ^-pulse stages of the interferometer sequence, with lattices acting as waveguides between these 
stages. The preceding analysis has developed the theoretical machinery for calculating phase shifts for these lattice 
interferometers. We now examine several of the most promising applications of lattice interferometers. 



A. Gravimetry and gravity gradiometry 

Lattice interferometers can be used to make extremely precise measurements of the local gravitational acceleration 
g. We proceed to calculate the phase shift for a lattice gravimeter. It is essential to note that whenever the two arms 
are addressed by different lattices they will be in different dressed state frames (where we recall that a dressed state 
frame is defined by the velocity of the corresponding lattice and by the distance that the lattice has traveled since the 
beginning of the interferometer sequence) . Let the velocities in the lab frame of the two lattices be denoted as (0 
and fLa™^(^)' respectively. The lattice velocities for the two arms in their respective dressed state frames will thus be 
^^Lauice(0 = «rab^(^) + 9^ ^ «o and w£™tico(*) = «Lab^(*) + 9^ " «o- We uote that vq is the velocity of the atom before 
the initial Bragg diffraction that splits the arms in momentum space, so that the two arms have different momenta 
after the Bragg diffraction. During the lattice loading period, the two lattices must be resonant (as described in Sec. 
Ill and in Appendix C) with the respective portions of the atomic wavefunction that they are addressing. Also, let 
Av{t) = WLab^(0 ~ '^Lab^(*) the vclocity difference between the two arms and Ad{t) = /J Av{t')dt' be the distance 
between the two arms. Where the interferometer sequence lasts for a time T, we can derive the phase shift for a 



lattice gravimeter using Eq. (13). We note that an additional contribution to the phase shift will arise if the two 



arms of the interferometer end up in slightly different dressed state frames. The arms begin in the same dressed 



state frame, and if Ad{T) = they will end up in the same dressed state frame, in which case Eq. ( 13 ) provides the 
final say on the phase shift. If Ad(T) differs slightly from zero, the two arms will end up in slightly different dressed 
state frames. Since we must compare phases in a common frame, it is convenient to boost both arms into the freely 
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falling frame using the transformation given in Eq. (A5) and the mathematical framework discussed in Appendix B. 
In addition to the phase difference from Eq. (13), we will then have an additional term in the phase shift equal to 
— jr^m{vf — I'o + gT)Ad{T), where mvf is the lab frame momentum at time T of a particular momentum eigenstate 
in the atomic wavepacket. After averaging, it follows from the discussion in [T that mvf will take on the value of 
the center of the momentum space wavepacket. We can combine the contribution from Eq. ( 13 ) and the contribution 
from boosting the two arms to a common frame to calculate the phase shift, where we also include a term 0Bragg to 
embody the net contribution to the phase shift arising from the Bragg pulses: 



Acj) 



Hvllt(t) + gt - vofdt - / -m«r (<) +9t- vofdt 



H<Titr~vllTmdt+ mgAv{t)tdt- mv^Av{t)dt 



l{Vf -Vo+ gT)Ad{T) + 0Bra 



- ^"^(^/ ^ "^'0 + gT)Ad{T) + 0Bragg 



(31) 



Integrating the second term by parts and simplifying yields: 



^^ = -nl, ^9^d{t)dt +-J^ ^m{vl^^\tr - vll^'itf)dt - jVfAd{T) + 0B.agg (32) 

When expressed in terms of lab frame quantities, it is apparent that Acj) contains terms corresponding to the propa- 
gation phase and the separation phase that typically appear in standard atom interferometer phase shift calculations. 
We note that in the dressed state frame, ideal Bragg pulses simply yield contributions in units of ±| to the overall 
phase shift between the arms, and these contributions can easily be made to cancel so that '/>Bragg = 0. However, we 
note that in some cases, corrections to the simplified picture of an ideal Bragg pulse due to such factors as gravity 
gradients, finite pulse and detuning effects (which can sometimes lead to a non-negligible propagation phase during 
the Bragg pulse), phase noise, or population loss may need to be considered. To avoid unnecessarily complicating 
our presentation, we will not present these corrections here. Instead, we emphasize that they are well-understood 
effects and refer the reader to other sources for further discussion pQ|5TJ[33]. Moreover, we note that these effects gain 
additional suppression for interferometers in a differential configuration so that they will often be below the mrad 
level P, as in the case of the gravity gradiometer discussed below. 

In the symmetric case where the velocities of the two arms in the lab frame are either opposite to each other or 
equal so that ^Lab^lO^ ~ ''^Lab^(^)^ ~ 0. Since we need the two arms of the interferometer to overlap at time T, it is 



convenient for us to choose Ad{T) = 0. Thus, the first term in Eq. (32 1 will constitute the only contribution to A0 



However, in an experiment, the parameters in Eq. (32) will undergo small fluctuations around their desired values 
from shot to shot, so that the other terms in Eq. (32) act as a source of noise. 

In order to cancel the effects of this noise, we can adopt a gradiometer setup in which an array of two or more 
gravimeters interacts with the same lattice beams. Although fluctuations in Ad{t) will still affect phase differences 
between gravimeters, which take the form — ^ /q rn{gi — g2)^d{t)dt, modern phase lock techniques will typically allow 
us to control the phase differences between the lattice beams well enough so that these effects are smaller than shot 
noise ^7\. When measuring a gravity gradient, the value of g will vary due to the gradient over the range of a single 
gravimeter. For linear gradients, we can calculate the phase shift in the presence of a gravity gradient by assuming 
that the value of g corresponding to the gravimeter is equal to its value at the center of mass position of the atom 
(see [3T] for a rigorous justification of this procedure) . When higher order derivatives of the gravitational field become 
sizable in comparison to the first derivative, this simple prescription may not suffice, and we can treat the problem 
perturbatively. If we want to measure an acceleration as well as a gravitational gradient in a noisy environment in 
which the fringes of the individual gravimeters are washed out, we can use dissimilar conjugate interferometers whose 
phase noise is strongly correlated as suggested in ^10^ . Appropriate statistical methods can then be used to extract 
the desired signal [38 . 

The effects on the phase shift of non-adiabatic corrections, lattice depth imbalances, and the finite spread of the 
atomic wavefunction in momentum space are considered in Sec. IV., Sec. V., and Appendix B. Based on the analysis 
in these sections, we conclude that a wide range of experimentally feasible gravimeter geometries exist that contain 
sufficiently adiabatic lattice depth and velocity ramps and that make use of symmetry in such a way that the net 
contribution of these corrections to the phase shift will be below the mrad level in a gradiometer configuration. We 
note that this paper has developed the mathematical machinery to calculate any such corrections to arbitrary precision 
if necessary. 
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When two lattices are used to manipulate the arms of an atom interferometer, one lattice will be on resonance with 
a given arm, while the other will be highly detuned. As long as we keep this detuning large enough and/or employ 
geometries with sufficient symmetry between the arms, the net effect of the off-resonant lattices on the final phase 
shift can often be made to be smaller than the mrad level in a differential configuration (e.g., a gradiometer) . For 
arm 1, the detuned lattice will manifest as an additional term in the discrete Hamiltonian, given by: 



detuned 



__6 „p-»0(t) 








_6 oe~*'^(*) 



_0 oe*'3(*) 





^ ^arm2'-^ 





V 





-fi op'''3(*) n 

(33) 

where (3{t) = - [a''™^(i') - a'^'"^'^{t')] Aur-dt' = -2kAd{t) and where f2arm2 is the lattice depth parameter corre- 
sponding to the detuned lattice that addresses arm 2 (32]. To obtain the correction to the Hamiltonian for arm 2, 
which comes from the detuned lattice addressing arm 1, we replace j3{t) with —l3{t) and riaim2 with riarmi- For large 
detunings, /3{t) will vary rapidly with time so that the contribution from the detuned Hamiltonian will be small due to 
the rotating wave approximation. Corrections arising from the detuned Hamiltonian can be solved for perturbatively 
using methods such as adiabatic perturbation theory. But we emphasize again that we can often avoid situations 
where this will be necessary. For example, we find from perturbation theory that to avoid population loss due to the 
off- resonant lattice as described in Sec. TIL, we should choose the off-resonant lattice to have a velocity that differs 
from that of the particular arm of the interferometer under consideration by an amount A?; 3> \ ft\vr (where is the 
depth parameter of the off-resonant lattice). We have verified this result with numerical simulations. However, in 
this regime, the off-resonant lattice can still sometimes cause a non- negligible energy shift, which we estimate with 
perturbation theory. The energy shift for arm 1 is AE'armi ~ f w~-^(?;r./Az;)^(riarm2)^- Analagously, the energy shift 
for arm 2 is Ai?arm2 ^ f w~^(z;i./Aw)^(riarmi)^- The relevant quantity in determining the correction to the phase 
difference between the arms is the difference in energy shifts Ai?armi — Aii'arm2, which is determined by the lattice 
depth imbalance between the arms. Where we let riarmi — ^ a-nd riaim2 = + Afi: 



h nAn , 1, /27r X 3.77kHz \ / 20 V f Vl V f AVl/n\ 

This result agrees with numerical simulations. Fig. [s] illustrates the dependence of (AE'armi — Ai?aim2)/^ on (Av/vr). 
For the same reasoning as in the discussion in the previous section, the net effect of lattice depth imbalances on 
the correction term treated here will often be further suppressed by orders of magnitude for interferometers in a 
differential configuration. Thus, as stated previously, in many cases the net effect of phase corrections due to the 
off-resonant lattices will be below the mrad level in a differential configuration. 

A lattice gravimeter can provide extraordinary levels of sensitivity. This sensitivity can be achieved over small 
distance scales by implementing a hold sequence in which the two arms are separated, manipulated into the same 
momentum eigenstate, held in place by a single lattice, and then recombined. In achieving compactness, we note that 
the fact that lattice interferometers are confined and can thus keep the atoms from falling under gravity during the 
separation and recombination stages of the interferometer as well as during the hold sequence is essential. Otherwise, 
for many configurations, the desired arm separation could not be reached without the atoms falling too great a distance, 
which would ruin the compactness of the interferometer. Hold times will be limited by spontaneous emission, which 
decreases contrast. Modern laser technology will allow us to use detunings of hundreds or even thousands of GHz, 
making hold times on the order of 10 s within reach [T5 . Gravimeter sensitivities using the hold method greatly exceed 
the sensitivities of light-pulse gravimeters while simultaneously allowing for a significantly smaller interrogation region. 
For example, for ^ 10^ atoms/shot and ^ 10^^ shots/s, a shot noise limited conventional light-pulse interferometer 
with a 10 m interrogation region can achieve a sensitivity of ~ 10~^^ g/Hz^/^. With similar experimental parameters, a 
shot noise limited lattice interferometer with a 10 s hold time and an interrogation region of 1 cm will have a sensitivity 
of ^ 10"^^ g/Hz^/^. If we expand the interrogation region to 1 m, we obtain a sensitivity of ^ 10"^'* g/Hz^/^. This 
remarkable sensitivity has a plethora of potential applications. Extremely precise gravimeters and gravity gradiometers 
can be constructed to perform tests of general relativity, make measurements relevant to geophysical studies, and build 
highly compact inertial sensors. Moreover, the fact that lattice interferometers can operate with such high sensitivities 
over small distance scales makes them prime candidates for exploring short distance gravity. One could set up an 
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array of lattice gravimeters to precisely map out gravitational fields over small spatial regions, as shown in Fig. |6] 
The knowledge obtained about the local gravitational field could be useful in searching for extra dimensions j39j as 
well as in studying the composition and structure of materials. 



B. Tests of atom charge neutrality 



Ultra-high precision gravitational measurements are certainly among the most promising applications of lattice 
interferometers, but the usefulness of lattice interferometers is certainly not limited to the study of gravity. By 
exposing the two arms of a lattice interferometer to different electrostatic potentials, tests of atom charge neutrality 
with unprecedented accuracy could be achieved |40j . The main advantage of a lattice interferometer in such a 
measurement is that the interrogation time can be significantly increased in comparison to the interrogation time 
achievable in a light-pulse geometry through the use of a hold sequence. Where thoid is the duration of the hold, V 
is the electrostatic potential difference between the two arms of the interferometer, e is the electron charge, and e is 
the ratio of the atomic charge to the electron charge, the phase shift is given by ^Vt^oid Based on the results of 
Sec. IV., Sec. V., and Appendix B, and assuming an identical configuration to that described in except for the 
inclusion of a hold sequence, we estimate that the phase error induced by the undesirable effects we consider will be 
below the proposed shot noise limit for this experiment (1 mrad). Any systematic phase error can be characterized 
by the methods we have developed. For a hold time of 10 s and for integration over 10^ shots, atom charges can be 
probed down to the region ofe^lO^^^. 
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FIG. 6; (color online) An array of lattice gravimeters such as that shown above can be used to achieve measurements of a local 
gravitational field with high spatial resolution. The trajectories shown hold in the lab frame. After separating the atoms in 
each gravimeter by a small amount, we can implement a hold sequence to greatly increase sensitivity. Bragg pulses are used at 
times to, ti, t2, and t^. Such an array could be used to study general relativistic effects, search for extra dimensions, examine 
local mass distributions, or measure Newton's constant. In addition, the ability of lattice interferometers with hold sequences 
to provide extremely precise measurements with a small interrogation region makes them ideal candidates for compact, mobile 
sensors. 



C. Measurements of — and of isotope mass ratios 

The ratio ^ is of particular interest because of its direct relation to the fine structure constant. Atom interferometry 
has previously been used to provide exquisite measurements of — . The most precise atom interferometric measurement 
to date was performed by Cadoret et at, who performed an elegant experiment combining Bloch oscillations and a 
Ramsey-Borde interferometer to measure the fine structure constant to within a relative uncertainty of 4.6 x 10~^ [llj . 



Eq. (32 1 indicates that if we apply different accelerations to the two arms of the interferometer, we will see a phase 
shift proportional to ^ that depends on the kinetic energy difference between the arms, which can be made extremely 
large. A differential configuration using conjugate interferometers (shown in Fig. [?]) could reduce the net contribution 
of such unwanted effects as laser phase noise and cancel the gravitational phase shift up to gradients ^U\. Such a 
geometry could provide an extremely precise measurement of ^ , as illustrated by the fact that we can achieve a phase 
shift of ~ 10^^ radians for a 5 m interrogation region and a 0.6 s interrogation time, corresponding to a shot noise 
limited sensitivity of ~ 10~^*^/Hz^/^ for the experimental parameters stated above. In this situation, the dominant 
unwanted effect would arise from non-adiabatic corrections, and the methods for calculating these corrections that 
are presented in the previous sections would need to be applied (we estimate a phase error ~ 10 rad). We note 
that even if the shot noise limit is not reached, the technique we have proposed could still improve over current ^ 
measurements. We emphasize again that to take full advantage of the sensitivity offered by lattice manipulations in 
atom interferometry, the methods we develop in this paper for calculating non-adiabatic corrections are absolutely 
essential. 
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FIG. 7: (color online) Conjugate interferometer geometry that could be used to measure ^. The trajectories shown hold in 
the lab frame. Bragg pulses are used at times to, ti, and t2 (at to, a sequence of multiple Bragg pulses is necessary to split the 
system into the two arms of the two conjugate interferometers). The phase shift of the lower interferometer is subtracted from 
the phase shift of the upper interferometer, which suppresses the effects of laser phase noise and eliminates the gravitational 
phase shift up to gradients [10]. With such a scheme, a measurement of to a part in 10^* may be possible, which could lead 
to the most accurate determination of the fine structure constant to date. 



The fact that all terms in Eq. (32 1 are proportional to m (except for the 4>BTa,gg term, which as we have explained, 
will often be negligible) can be exploited to provide high-precision measurements of isotope mass ratios by using an 
interferometer geometry in which the two isotopes follow identical trajectories. Isotope mass ratios could be relevant 
to studies of advanced models of the structure of the nucleus [H] . Conversely, if we know the mass ratio of two isotopes 
sufficiently well, we can use such a geometry to precisely measure accelerations, where the two isotopes provide two 
dissimilar conjugate interferometers. The phase noise of these two interferometers will be extremely well correlated 
because they are topologically identical. This also eliminates the need for the additional lattice beams required to 
form the second, topologically distinct interferometer that would be needed if we only had a single isotope. Note 
that this scheme to construct dissimilar, topologically identical conjugate accelerometers would not be possible for a 
light-pulse geometry, for the leading order phase shift of light-pulse accelerometers is independent of isotope mass. 



D. Gyroscopes 

Lattice interferometers can also be used to build compact and highly sensitive gyroscopes. There are multiple 
possible schemes in which optical lattices can enhance gyroscope sensitivity. One such scheme is to modify a typical 
atom-based gyroscope by replacing the Raman pulses with LMT lattice beam splitters, increasing the enclosed area 
of the interferometer and hence its sensitivity to rotations. The phase shift can be written in Sagnac form as 
■ A, where (l is the rotation rate vector and A is the normal vector corresponding to the enclosed area of the 
interferometer [1]. The gyroscope described in [6] achieves a sensitivity of 2 x 10^^^/Hz^/^. Replacing the Raman 
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pulses in this experiment with 200fi.fc lattice beam splitters would increase the sensitivity by a factor of 100. In 
such a configuration, we estimate that each beam splitter could introduce a non-adiabatic phase error of ~ 1 rad 
if the arms of the interferometer are not split symmetrically. However, beam splitter configurations that exploit 
symmetry between the arms can reduce this effect by orders of magnitude. Another option is to use optical lattices 
along multiple axes to provide complete control of the motion of the atoms in two or three dimensions (this control 
is only achieved in the region in which the lattices overlap, necessitating the use of wide beams). Analagous to a 
fiber-optic gyroscope, the atoms could be guided in repeated loop patterns, with the two arms rotating in opposite 
directions. Geometries in which atomic motion is controlled in multiple dimensions could also expand the possibilities 
for other applications of lattice interferometry (such as measurements of gravity) by allowing for the measurement of 
potential energy differences between arbitrary paths. For instance, a compact array of three orthogonal lattice gravity 
gradiometers could be used to measure the nonzero divergence of the gravitational field in free space predicted by 
general relativity [H]. 

VII. CONCLUSION 

We have presented a detailed analytical description of the interaction between an atom and an optical lattice, using 
the adiabatic approximation as a starting point and then proceeding to rigorously develop a method to calculate 
arbitrarily small corrections to this approximation using perturbative adiabatic expansion. We have applied this 
theoretical framework to calculate the phase accumulated during a lattice acceleration in an LMT beam splitter. And 
we have proposed atom interferometer geometries that use optical lattices as waveguides and discussed applications 
of such geometries, using our theoretical methods to add rigor to this discussion. We are working toward the ex- 
perimental implementation of lattice interferometers and LMT lattice beam splitters, and we hope to explore the 
applications we have discussed. In this experimental work, we realize that we will have to contend with a number of 
unwanted systematic effects, such as spatially varying magnetic fields, imperfections in the lattice beam wavefronts, 
and inhomogeneity of the lattice depth across the atomic cloud. We have studied these effects using both analytical 
and numerical methods, and we are optimistic that they can be significantly mitigated for a wide range of experi- 
mental parameters-a conclusion that we hope to verify experimentally. Many of the unwanted systematic effects that 
are relevant to lattice interferometers are also shared by light-pulse interferometers and can therefore be dealt with 
using similar methods. Therefore, we believe that it will likely be possible to realize lattice interferometers in existing 
apparatuses originally constructed with light-pulse geometries in mind. 
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APPENDIX A: BOOSTING BETWEEN DIFFERENT FRAMES 

For the purposes of this paper, we consider unitary transformations that consist of a translation in position space, 
a boost in momentum space, and a time-dependent change of phase. Such a transformation has the general form: 

U{t) = eS'*WPe~^™^(*)*e^''(*) (Al) 

We note we are free to choose the boost parameters d{t) and v{t) arbitrarily. That is, we do not have to choose them 
so that v{t) is the rate of change of d{t). The operators x and p transform under U{t) as follows: 

(l{t)xU\t) = x + d{t) 
U{t)pU\t) = p + mv{t) 

(A2) 

We now consider the Hamiltonian in a frame that is freely falling with gravity, which takes the form: 
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p2 

2m 



kx - k{Di^,b{t) + ^gf) 



(A3) 



We can transform from the freely falling frame to the lab frame by applying the appropriate Galilean transfor- 
mation UpFit), which corresponds to specifying d{t) = hgt^, v{t) = gt, and 6{t) = ^mg^t^, so that HLab = 



Since the only position dependence in H^y comes from the lattice potential, we could readily approach the problem 
of finding the dynamics of the system by working in the freely falling frame. However, it will prove to be useful from a 
calculational standpoint to transform to a third frame, with a Hamiltonian resembling that describing the atom-light 
interaction from the point of view of dressed states. We note that although we could have performed a boost directly 
from the lab frame to the dressed state frame, it is useful to introduce the freely falling frame for pedagogical reasons, 
since it is the frame in which calculations for atom interferometry are typically performed. In Appendix B, we use 
the transformation between the freely falling frame and the dressed state frame to highlight the parallels between a 
lattice beam splitter and a typical light pulse beam splitter. We absorb the initial velocity uq of the atom in the lab 
frame (and hence also in the freely falling frame) into the dressed state frame, so that velocity vq in the lab frame 
corresponds to velocity zero in the dressed state frame. The Hamiltonian in the dressed state frame is: 

Fds = - («Lab(i) +9t- va)p + 2hn{t) sin^ (kx) (A4) 
2m 

The unitary transformation Uds that transforms from the dressed state frame to the freely falling frame, so that 
Hff = UDsit)HDsUj)g{t) + ih (^§iUDs{t)^ ^Ds{t)^ corresponds to boost parameters d{t) — — (I?Lab(i) + ^5^^) and 
v[t) = —vq with a time-dependent phase factor: 



U-Ds = es™"''*e""s^-°^°''**)+5^*^^^e""s(5™''o*) (A5) 

The ability to boost to a frame in which the Hamiltonian contains no position dependent terms outside of the lattice 
potential is contingent upon the assumption that the external potential in the lab frame (not including the lattice 
potential) is linear in x. However, real-world potentials such as the potential corresponding to Earth's gravitational 
field will deviate somewhat from this assumption. Any such deviations would manifest as residual position dependent 
terms in the dressed state Hamiltonian, which we collectively refer to as V . Under the semiclassical approximation, 
we neglect the effects of V on the time evolution of the atomic wavepacket. This approximation is valid when the 
energy scale of V over the spread of the atom's wavefunction (which is on the order of magnitude of the expectation 
value of V in the atomic wavepacket) is much smaller than the energy scale of the lattice potential and is small relative 
to the time scale of the experiment, which is the case for a wide class of experimental parameters. For example, in 
the case of a Rubidium atom wavepacket with a spatial spread of cr ~ 100/xm in the gravitational field at the Earth's 
surface (which has a gradient of A ^ 10^^ s"^), the energy scale of V will be ^ imAcr^ ~ /i x (10^® Hz). This energy 

scale is smaller than that of a lattice of typical experimental depth (^ SE",-, where Er is the recoil energy %^) by a 
factor of ^ 10^° and is small on a time scale of ~ 10 s. The effects of linear gradients and of more general potentials 
can be accounted for through a straightforward generalization of the results presented in this paper, as discussed in 
greater depth in [3T] . 



APPENDIX B: GENERALIZING TO THE CASE OF A FINITE WAVEPACKET 

In Sec. II., we discretized the Hamiltonian using the basis of momentum states \2nhk) for integer n. We now 
generalize our results to the case of a finite wavepacket. 

Throughout our analysis, we have worked mainly in the dressed state frame, since this frame is particularly con- 
venient for describing phase evolution in a lattice. For interferometers where we use lattices as waveguides for the 
atoms, we will want to perform the entire phase shift calculation in this frame. However, phase shift calculations for 
light-pulse atom interferometers are often performed in the freely falling frame. Thus, for applications where lattice 
manipulations are used for beam splitters and mirrors in a light-pulse geometry, it is useful to convert the phase 
evolution we calculate in the dressed state frame to the freely falling frame. We thus present the general results 
derived in this appendix in the freely falling frame. 
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In Sec. II., we considered a particular dressed frame in which the initial velocity of the atom is boosted to zero. We 
now introduce an unboosted dressed state frame that is related to the freely falling frame by a translation in position 
space with no boost in momentum space, so that the unitary transformation Uus,-, — e^'s(-^i''»'=(*)+29* )p transforms 
from the unboosted dressed state frame to the freely falling frame. The Hamiltonian in this frame is: 

^DSo (wLab(i) + gt)p + 2nn{t) sin^ {kx) (Bl) 

Now, say that before the lattice acceleration, the state that we are accelerating is described by I^'pfC^o)) in the freely 
falling frame. We can then transform this state vector to the unboosted dressed state frame, describe its evolution to 
the final time in this frame, and transform back to the freely falling frame. Where Tdso(^':^) is the time evolution 
operator the takes us from time t to time t' in the unboosted dressed state frame, we can write: 

|*FF(i/)) = %So(i/)Tbs„(i/,io)C>iso(*o) l*FF(io)) (B2) 
Denoting the initial momentum space wavefunction in the freely falling frame as ^ff(p, ^o) = (pI^ffC^o)); we can 



express Eq. (B2) as 



\^Mtf)) = Umoitf)fBSo(.tf,toMs„{to) j dp\p)^MpM 

dpJ7DSo(i/)T'DSo(i/,io)b)e*(^^^^(*«'+^^*«)P'fFF(p,to) 

(B3) 

where we have used the linearity of the operators to bring them inside the integral. We now consider how each 
momentum eigenstate \p) evolves in the unboosted dressed state frame. That is, we must calculate Tbso(*/j^o) \p) 
for each p. In order to do so, we introduce a class of boosted dressed state frames DSp parameterized by p, so that 
momentum p in the unboosted dressed state frame (which is just the frame DSq) corresponds to momentum zero in 
the frame DSp. In essence, where u = ^, frame DSp travels with velocity v with respect to the unboosted dressed 
state frame. In frame DSp, the Hamiltonian takes the form: 

^^DS, = (i^Lab(t) +gt- v)p + 2hn{t) sin2 {kx) (B4) 

^ 2m 

where we note that the unitary transformation that transforms from frame DSp to the unboosted dressed state frame 
is: 

Up{t) = e^P*e~*[~*-°^"'='^*^+5^*^^™'""''5™^'^*] (B5) 
Where TjjSp{tf ,^0) is the time evolution operator in frame DSp, we can write: 

^DS„(t/,to)b) = JJp{tf)fus,itf,to)Ulito)\p) 

(B6) 

The problem of determining the evolution of the momentum eigenstate \p) in the unboosted dressed state frame 
thus reduces to evolving the momentum eigenstate |0) in the frame DSp, which we know how to do from the preceding 
sections. We assume that the momentum space wavefunction is narrow enough so that the momentum eigenstates 
we consider are resonant with the lattice (that is, as we recall from Sec. III., the magnitudes of their velocities with 
respect to the lattice are less than Vr). Where the lattice acceleration is chosen so as to transfer a momentum of 2nhk 



to the atom, the result from Sec. III. tells us that the time evolution in Eq. (B6) yields: 



Tbs,(i/,io)|0) =e'^-|2nfifc) (B7) 
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for phase: 



r,^ , <atticc(*) 

^" -'to 



(B8) 



where v^^^.^.^^^{t) = t;Lab(^) + gt — v is the lattice velocity in frame DSp. For the sake of pedagogy, at the moraent 
we neglect the phase arising from the lattice depth, the phase arising from the small periodic variations in the 
lowest eigenvalue, and the phase arising from non-adiabatic corrections, where we note that we could calculate these 
contributions if needed. Evaluating the integral in Eq. (B8), we obtain: 



2h 



(B9) 



Substituting this result into Eq. (B6), applying the transformation Up{tf), and canceling terms in the exponentials, 
we find that: 



(BIO) 



We can now evaluate Eq. (B3), which gives us the final result: 



(Bll) 



where: 



b) = ^ / (fLab(0 + gtfdt - 2nk{Di^,^{tf) + -gtj) - |AD (B12) 

J to 



ypF 



for AD = (-DLab(t/) + yt}) - {Di^^M + ^94)- 



Observe that the second term in Eq. (B12) is what would typically be called the laser phase in a light-pulse atom 
interferometer for an n-photon beam splitter. We note that as long as the resonance condition is met for the momentum 
eigenstates we are considering so that they are accelerated, the phase evolved during a lattice acceleration in the DSq 



dressed state frame, is independent of p, as shown in Eq. ( B10[ ). (This is true up to non-adiabatic corrections and the 



small periodic variations in the ground state eigenvalue) . This makes dressed state frames particularly convenient for 
performing calculations involving wavepackets. In contrast, in the freely falling frame, the accumulated phase (jj^F is 
dependent on p, but this dependence cancels in the final expression for the phase shift between the two arms of an 
interferometer as long as the distance travelled by the atom while locked into a lattice is the same for both arms. 
Note that momentum dependent contributions to the total phase shift must arise when we treat the problem purely 
in terms of dressed states, since the total phase shift is an observable quantity and must therefore be independent of 
the frame in which it is calculated. The key point to realize is that if the total distance traveled in the lattice is not 
the same for both arms, then the two arms will end up in two different dressed state frames. 

We now consider the effect on phase evolution of the periodic term po{a, fl) in the expression for the ground state 



eigenvalue of the dressed state Hamiltonian given in Eq. (11). This term leads to a momentum dependent correction 



f)p to the evolved phase 0p described in Eq. (B8|, where: 



^ /" „ ^Lattice (^) rS/ 



to 



We note that po ^ "Latticc(*) ^ ^(^)^ can be calculated using the truncated matrix approximation described in Sec. IV. 

Since po ^^^i^siiiss^ ^ ^{t)^ is periodic in i'Lattico(0 with period 2vr, the contribution S(j)p™^ — (50^'™^ of this correction 
term to the phase difference between the two arms of an interferometer will be highly suppressed if both arms load 
the atom into the lattice near the center of the zeroth band (as discussed in Sec. III.) and undergo a nearly integral 
number of Bloch oscillations so that the effects of this periodic variation in the ground state eigenvalue will be largely 



common to the two arms, as verified by estimation of the integral in Eq. (B13) and by numerically solving the 
Schrodinger equation. 
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APPENDIX C: THE RESONANCE CONDITION FOR AN ATOM TO BE LOADED INTO THE 

GROUND STATE OF A LATTICE 



Here, we derive the resonance condition stated in Sec. III., which says that for a sufficiently slow lattice depth ramp, 
the atom will be loaded into the ground state of the lattice if the initial velocity of the lattice is within Vr of the initial 
velocity of the atom. First, we examine the eigenvalues of H when ^{t) — 0. For Cl{t) = 0, is diagonal matrix. Since 
H is a, matrix defined in terms of the basis of momentum eigenstates \2nhk) for integer n, when H is diagonal these 
momentum eigenstates are also the eigenstates of H , with corresponding eigenvalues £"„ — 4E,. (n'^ — rta(i)) (which 
are just the diagonal elements of H). We note that the effect of the — (wLab(^) + gt — vo)p term in the Hamiltonian is 
to encapsulate the dependence of the eigenvalues and eigenvectors of the Hamiltonian on the lattice velocity for any 
given lattice depth, which can be best conceptually understood from the point of view of the dressed state picture 
explained in |30]- The resonance condition states that if the initial state is \2nQhk), then the initial velocity of the 
lattice must be within Vr of 2noVr- Thus, the initial value ao of the dimensionless velocity a must be in the range 
(2no — 1, 2no + 1). We will now show why this condition on ao implies that the ground state of H with 17 (t) =0 is 
\2nQhk). Where we write ao = 2no + S and let n = uq + An for integer An, the eigenvalues of H with the lattice 
depth set to zero are: 



{no + Anf - {no + An) {2no + 5) = 4£^r [-'^o + (An)^ - noS - An5] 



(CI) 



For 1^1 < 1, En will be most negative when An — 0, so \2nohk) will indeed be the ground state of the Hamiltonian 
when n{t) = 0. Therefore, as long as we ramp up the lattice slowly enough and keep the lattice velocity constant 
while doing so, the adiabatic theorem tells us that the atom will remain in the lattice ground state throughout the 
ramp up process provided that the ground state never passes through a point of degeneracy, which is indeed the case 
for fixed a and \S\ < 1. 

The maximum rate at which we can increase the lattice depth while still maintaining the validity of the adiabatic 



'MH\0) 



(£l-£o)^ 



where 



approximation depends on |i5|. This statement follows from the adiabatic condition 

|0) and |1) respectively denote the ground state and the first excited state of the Hamiltonian. When the resonance 
condition holds, note that the first excited state will be |2(no + sign{S))hk) . As \6\ increases toward 1, the gap between 
En„ and Eng+sign(S) (i-e. the energy gap between the zeroth and first bands) decreases. To see why this is true, we 
observe that: 



i?«o+sig„(5) - = AEr{l - sign{S)S) = 4Er{l - \5\) (C2) 

Therefore, for the adiabatic approximation to hold, the maximum rate at which we can ramp up the lattice decreases 
as we increase \6\. This result is identical to the statement in Sec. HI. that adiabatic loading is more difficult to 
achieve near the border of the first Brillouin zone. 

After ramping up the lattice, we accelerate it until its velocity is within Vr of 2nfVr where \2nfhk) is the target 
momentum eigenstate. The ground state of the Hamiltonian when the lattice has finished ramping down will therefore 
be \2nftik), and as long as the entire process is adiabatic so that the atom remains in the ground state, the atom 
will indeed end up in the target state. By analogy to the previous discussion, the maximum rate at which we can 
adiabatically ramp down the lattice increases as the difference between the final lattice velocity and 2nfVr goes to 
zero. 



APPENDIX D: PERTURB ATIVE ADIABATIC EXPANSION AT HIGHER ORDERS 

In order to calculate non-adiabatic corrections at arbitrary order, we define Ca^t^ >m^n^j recursively in the 

natural way based on the notation of Sec. IV: 

Ca^b^...^„r^n^j ^ - ^^b^. . (i') ^^^e^'^" (O"'., (*')] rf^' (Ql) 

For arbitrary order fc, we can then write: 

h^j'\t) = Cj + ^ Cn^j + ^ ^ Crn^n^j + ' ' ' + ^ ^ ' ' ' ^ Ca^b^ >m^n^j (D2) 
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where the last term includes k sums and Cj = bj . In order for the expansion to be practical from a calculational 
standpoint, the series must converge quickly enough so that we do not have to find an unreasonable amount of terms. 
There are two types of terms that we can often neglect. First, where we assume that the system is initially in the 
ground state, it is often possible to neglect many terms of the form Co-i-n so that we only need to examine a relatively 
small number of eigenvectors. Terms of this form indeed decrease rapidly as n increases because the overlap of the 
nth eigenstate with the zcroth eigenstate is essentially nonexistent for large enough n. Second, we can often neglect 
all terms of order greater than some cutoff value. We will be able to do so as long as the time scale over which we 
are solving the problem is small enough so that integrating a fcth order correction term against a factor of the form 

yields a {k + l)th order correction term that is much smaller than the correction term of order k. 
Note an important nuance in how we have phrased the above condition~we have mentioned nothing about a 
Hamiltonian that varies slowly in time. The convergence of the perturbative series depends on the time interval of the 
solution. As long as the Hamiltonian does not vary at an infinitely fast rate, we can always work on a small enough 
time scale so that the series converges rapidly. To solve the problem on time scales for which the series does not 
converge quickly, we can simply break the problem into multiple parts. This method provides with us with a means 
to describe the system for a Hamiltonian that changes arbitrarily fast in time. Having a slowly varying Hamiltonian 
just serves to allow us to solve the problem without dividing it into as many parts (much of the time we will not have 
to divide the problem at all, and in Sec. V. we derive conditions for when this will be the case), thus making the 
calculation significantly easier. 
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